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1. Motivation 



Since the early days of quantum field theories, it has been argued that their perturbative weak 
coupling expansions are asymptotic [pp. Usually, the asymptotic behaviour is identified using semi- 
classical methods such as instantons. Asymptotically free theories with marginal operators consti- 
tute a special case. Examples of such theories are four-dimensional non-Abelian gauge theories or 
the two-dimensional 0(N) model. For these the structure of the operator product expansion (OPE) 
is believed to give rise to a specific pattern of asymptotic divergence known as (infrared) renor- 
malon [Q, |3j] that cannot be obtained using standard semiclassical methods. In the two-dimensional 
0{N) model renormalon effects were found (though suppressed) in an explicit calculation in the 
large N expansion [Q]. For four-dimensional non-Abelian gauge theories no such proof exists. On 
the contrary, the possibility that the renormalon either does not exist or is very small has been raised 
in recent years, see e.g. [M 01. 

At present, the only way to unambiguously settle this issue is by performing weak coupling 
expansions to sufficiently high orders in the strong coupling parameter a. In particular, in SU(3) 
gauge theory, the problem has been addressed by simulations in lattice regularization, using Nu- 
merical Stochastic Perturbation Theory (NSPT) [0-0- These simulations, that mainly focused on 
computing the plaquette and small Wilson loops to perturbative orders as high as a 20 [|To|], so far 
had no success in uncovering a renormalon. 

The situation radically changed recently when, for the first time, accurate agreement between 
the explicit computation of perturbative coefficients and the expectation from the renormalon anal- 
ysis was reported [11]. Here we discuss some aspects of this analysis where the comparison was 
performed using another observable, the static self-energy. This, on theoretical grounds, is far more 
suited for a renormalon study. The simulations were conducted on a large set of different lattice 
volumes, for orders as high as a , and finite size effects were carefully investigated. 



2. Renormalon Theory 

In QCD, the expansion of a generic observable K as a power series in the coupling a, 

K = 1 £k n a n , (2.1) 

n 

is believed not to be convergent but at best to be asymptotic. Technically, renormalons appear as 
singularities in the Borel plane, making a Borel summation impossible. The way in which the 
series K diverges is closely tied to the OPE, commonly giving rise to terms k„ ~ a n d n\, where ad is 
a constant. For small orders n, the successive contributions k n a n to the series reduce in size down 
to a minimum at no l/(\a c i\a). The series should be truncated at no and inevitably one has to 
deal with an ambiguity of the order of the minimum term, k nQ a n ° ~ exp[— l/(|a^|a)]. 

The OPE allows one to isolate the short (given by Wilson coefficients d(q,n), where i is 
the dimension) and long distance effects (matrix elements (0,-(/i,A)}) of an observable R(q,A) by 
means of a factorization 

R = C (q,fl)(O (fl,A))+C d (q,fl)(O d (fl,A))(^\ +•••. (2.2) 
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Here, /I is the factorization scale separating the perturbative and low momentum scales q and A 
from each other: q 3> }i ^> A. 

For the plaquette, (Oo) = 1, followed by the dimension d = 4 gluon condensate as the next 
higher dimensional non-vanishing operator. The lower bound for the accuracy of expanding Co is 
of 0(h A /q A ) ~ k nQ a"°, the size of the minimum term, as can be seen from 



exp 



1 



where \a&\ 



ft = n 



(2.3) 



\ad\cc J 2nd 
The (infrared) renormalon ambiguity of this perturbative series cancels that of the next order non- 



perturbative matrix element in eq. fl2.2| ). Hence, the physical quantity R is well-defined. It follows 
that series expansions with the smallest d (and therefore n ) are those in which one should be able 
to detect renormalons most easily (for d = 1 four times "faster" than in the plaquette case). Among 
such candidates are the pole mass or the associated self-energy of a static source, 



8m = -Y J c n a n+l {\/a). 



(2.4) 



n>0 



a 1 , the inverse lattice spacing, serves as a UV-regulator. The large n behaviour of the coefficients 
c„ reads 



„_>. , p y r(n + i+b) ( b 



(2.5) 



where the coefficients b and s\ can be found in Jl2| ] and N m is a normalization constant that cancels 
when constructing ratios 



c„ 1 



c„-i n 



271 



1 + - 



(2.6) 



3. Numerical Simulation 

The simulations are performed within the framework of NSPT [^-01 that allows for a direct 
determination of perturbative series coefficients. A key ingredient of NSPT is the Langevin update, 
for which we employ the ff(Ar 2 ) integrator described in [13, 14]. Its convergence behaviour 



regarding the time step At is remarkably flat, so that we can circumvent a Ar-extrapolation and 
work at a fixed value At = 0.05. We implement periodic boundary conditions in time and twisted 
boundary conditions [ 15 - [l8|] in the three spatial directions, completely eliminating zero modes. 
Twists in two directions would have achieved this, too [|I^], but with reduced numerical stability. 
Employing the Wilson gauge action, we calculate the temporal Polyakov line 



lW(N s ,N t ) 



1 



i 

•~d R ' 



-tr 



N T -l 



n f?w 



«4=0 



(3.1) 



on hypercubic volumes with Ns and Nj lattice points in spatial and temporal directions, respec- 
tively. R denotes the representation of the link U^(n), of dimension d^, and we study both, the 
triplet and the octet cases. We repeat the measurements using stout-smeared [ p0| ] (smearing param- 
eter p = 1 /6) links in the temporal direction. Altogether, this amounts to four distinct self-energies 
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Figure 1: Comparison between the global fit and data for n = 9 {6"{a )) 



(allowing for checks of Casimir scaling and universality) whose perturbative coefficients c n are 
linked to the Polyakov line via 



P(N S ,N T ) = - — > dm. 

aNr 



(3-2) 



The dependence on Nt and Ns can be parameterized flllj] as 



aP = £ \c n a" +l (a- 1 ) - ^a^aNs)- 1 ) + ■ 



n>0 



A (1) 



A (2) 



£ c n +Al 1) (N s ) + K 2) (N s ,N T ) a n+ \a~ x ), 

1 



[/ n + logs{(AT s )] , 



N s 

{ Vn ~ W s [V " + 8Vn + logS " ( ^ )] } + N 1 { Wn ~N~ s [Wn + 8Wn + l0gS " W] 



(3.3) 



For sufficiently large Nj, the correction A„ dominates over A» , which comprises the leading 
/Nj, 1 /Ng) lattice artefacts. In contrast to the l/Nj- correction terms, the l/N§ corrections do 
not have a significant effect on our fits. Therefore, we neglect the terms containing w n and Sw„. 



The A„ contribution arises from interactions with mirror charges on lattice replica [BTJ, which 
result in an effective static potential between charges that are separated by distances aNs, but with- 
out self-energies. Consequently, the high order behaviour of the coefficients c n and /„ depends on 
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Figure 2: The ratios c„/(«c„_i) for the smeared and unsmeared fundamental static self-energies, compared 
to the prediction eq. at different orders of the l/n expansion. 



the very same renormalon. Precisely for this reason one needs to disentangle dm from the 1 / (aNs) 
correction carefully, in order to isolate the d = 1 renormalon. 

The log-terms in A„ are due to the running of a according to the j8 -function. In the lat- 



tice scheme the first three /3 -function coefficients /3o, j3i and /?2 have been calculated [22]. The 
lack of knowledge of higher j8y is leviated by the fact that these do not multiply the leading 
fn,fn-l)fn-2,fn-3 coefficients: at sufficiently high n, the /„ grow factorially and therefore sub- 
leading terms become strongly suppressed. 

Eq. ( j3j| ) is used for a global fit to all orders and all data from all our lattice volumes (see Ta- 
ble I of [|IT[]). To illustrate the adequacy of the ansatz, in fig. |] we compare the fit to the unsmeared 
triplet data, for the case of n = 9. Note that the ansatz requires only four parameters per order, 
yet it results in reasonable values # 2 /jVdf ~ 1-29 and 1.46 for smeared and unsmeared triplet data, 
respectively. 

Our main finding is shown in fig. ||, where we confront the infinite volume extrapolated data 
for c n /(nc n -\) with the theoretical prediction, eq. (2.6). Since smearing is a local operation that 



has little effect on long-range physics, the smeared and unsmeared data should exhibit the same 
large n behaviour, that is dictated by the infrared renormalon. This universality is supported by our 
data, as fig. illustrates. The octet representation data are also fully compatible with these findings, 



as we will detail in [14]. 
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Figure 3: Contributions to the sum eq. (2.4) for three different values of a — 1/ (|flrf|no), chosen such that 
no =7,9 and 12. 



In order to visualize the asymptotic behaviour of the perturbative expansion, in fig. |3] we dis- 
play the order-by-order contributions to the sum eq. (2A). For no = 1,9 and 12, three different 
values of a = 1 /(|«d|«o) are chosen (note that «o = 7, corresponding to j3 = 3/(27ia) ps 5.85 or 
aT x ps 1.6 GeV, is within the range covered by non-perturbative lattice simulations). The contribu- 
tions indeed become minimal at orders close to the predicted values no ps n + 1. 

With the infinite volume coefficients c n at hand, we can compute the normalization of the 
pole mass renormalon, see eq. (2.5). For smeared triplet data we get Af'f = 18.6(4), compared to 
N)n = 19.0(3) for the unsmeared case. As A^N^ S = A M Nj,f, this translates to 0.65(2) in the MS 



scheme. This agrees remarkably well with the estimate, N^ s ps 0.62 [123J, 124], from an expansion 



in the MS scheme up to &{a?). Hence, the d = 1 renormalon manifests itself in a lower bound of 
~ 0.65 Ajyjg for the precision with which one can determine the heavy quark pole mass. 



4. Summary 

Within the framework of NSPT, we have computed the static self-energy of SU(3) gauge 
theory in four spacetime dimensions to ^(a 20 ) in the lattice scheme. Simulations on various lattice 
volumes and a careful treatment of finite size effects are vital for our main finding, the onset of a 
factorial growth of the self-energy coefficients at an order n ps 9, in agreement with the theoretical 
renormalon prediction. 
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